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ABSTRACT 


An automatic-throttle compensation system has the capability 
of removing the problem of aircraft speed instability in the carrier 
landing approach. Various systems have been proposed and tested, each 
differing in input sensed variables. Two representative systems are in~ 
vestigated analytically and by the use of digital programmed Nyquist 
plots and analog simulations. An attempt is made to determine optimum 
gain constants by using various forcing functions in the analog simulation. 
Comparisons of the responses of the two systems are made using time 
history analog records. A digital computer stability program, applicable 
to any aircraft, is included. 
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AN INVESTIGATION OF THE EFFECT OF AUTO-THROTTLE DEVICES 

ON AIRCRAFT CONTROL IN THE CARRIER LANDING APPROACH 
i Introduction 

The carrier landing approach speed of current swept wing aircraft 
is influenced by ee factors. One of the most important of these factors 
is often referred to as speed instability. Speed instability can be ex- 
plained with the help of Pig. 1. It will be noted that for Aircraft "A", a 
slight deviation in airspeed from the normal approach speed makes very 
little difference in thrust required. This curve is typical of the AF-l1E 
type aircraft. On the other hand, a deviation towards slower speed for 
Aircraft "B" causes a large increase in thrust required. This increase, 
if uncorrected, causes a loss of altitude. Thus, in effect, the normal 
approach of Aircraft "B" is made on the "backside" of the thrust versus 
airspeed curve. This curve is typical of the F-8 aircraft. Essentially 
this steep backside of the curve is caused by a high value of Cam ; 
drag increase due to angle of attack (or lift). 

Still another factor contributing to airspeed instability is slow 
engine response. With a comparatively long time lag in engine response, 
any hesitation on the part of the pilot in correcting an incipient airspeed 
loss leads to wide variations in altitude or airspeed before effective 
thrust can be realized. 

With an aircraft that displays speed instability tendencies, either 
a faster approach speed or an automatic compensating device is required 


to prevent undesired altitude loss due to airspeed deviations. Faster 





approach speeds are not desirable for carrier approaches. It will be the 
purpose of this paper to investigate some methods of applying automatic 
power compensation to counteract airspeed instability. A swept wing 
fighter,, patterned after the F-8 airplane is used as the vehicle of the 
investigation. Both digital and analog computers are exploited in simu- 
lating auto-throttle and anti-drag system response to various inputs. 
Real-time, pilot controlled, analog simulated mirror approaches are also 
used in qualitatively evaluating each system. 
This investigation is the first of a continuing series to be conducted 

at the U. S. Naval Postgraduate School to evaluate systems of automatic 


airplane control. 





ae Discussion 

The aircraft-carrier mirror-landing approach is made at a constant 
speed and a constant angle of glide slope. This type of approach allows 
precision control of the landing, a factor of the utmost importance in 
carrier operations. However, the approach is usually made at a minimum 
airspeed that is dictated by performance and/or aircraft flying qualities. 
Thus, the airspeed must be monitored closely while making flight path 
angle adjustments. A simplified equation (given in Ref. 1) shows that 


flight path angle, YM , 1S approximated by the formula: 


yer~l+dVd ca’ a 


Ee. (1) Ww, aa 9 


Flight path angle, then, is dependent on lift/drag (L/D) ratio, thrust/ 
weight ratio, and rate of speed change. In the mirror approach the pilot 
can adjust ¥ and airspeed by manipulating the throttle and elevator. The 
manner in which each of these controls can be used, however, is depend- 
ent upon the L/D ratio. For aircraft with a high L/D ratio, the elevator is 
used to control 6 , While thrust is used to maintain constant airspeed. 
With a L/D ratio in the vicinity of 4 - 5 the rate of speed change is large 
during maneuvers. If, in addition, the L/D ratio decreases with increasing 
Cr in the approach speed range then speed instability is accentuated. 

A variation of lift/drag ratio versus lift coefficient for the F-8 air- 
plane is shown as Fig. 2. It can be seen that at the approach speeds 
used, (1.1-1.2 Voz ), L/D decreases rapidly with increase in Cr: 


As reported in Ref. 1, pilots stated that in order to change & for an 








airplane of this class, increased reliance on thrust rather than elevator 
was necessary. In the case of the F-8 airplane, however, height response 
to throttle input is slow due to small thrust line angle of attack and 
negligible trim changes due to power. The F-8 pilot, therefore, is forced 
into using the stick as the primary flight path control. Then, since the 
airspeed mode of oscillation (the phugoid) is lightly damped , the aircraft 
will hunt about a new equilibrium airspeed and glide slope angle in res- 
ponse to an elevator deflection. These speed changes in the 5 to 10 
seconds following elevator input are disturbing to the pilot. 

Another way of interpreting airspeed instability is referred to in 
Ref. 2. A curve similar to Fig. 1 can be drawn for thrust required versus 
dynamic pressure, q. A definition of a stable speed regime can then be 


postulated as that portion of the curve where 


ol Treguired 
am 


To express this inequality in terms of the lift and drag coefficients, 


Ineq. (2) 


Foe 


we use the fundamental equilibrium equations for level flight: 
Cy Sq = Drag = Treg 
en S¢. = Lift = W~ 


therefore li ayy Cp’ 
re = 
or 
ane 
and a Cy 


and A leg ais (Ww  ) = w al ( =) 








since W and S$ are constant. 





Therefore oA Treg = 3 d ( 32 sd() 
ol $ Tey 


c é 
which result is substituted into Ineq. (2). S ol ( S| 
ete 
A Cu ) 
Since 5 is greater than O, this implies that 
al ( Co ; 
Cy dl (onl 7° 
Taking the differentials of numerator and denominator: 


= Cy dlr * eft fb - ¢,dlo' - Cr’ dl Cr 
“ex ACy 7 -~ACy 


cL 


?O 


7 oO 








or 


Ineq. (3)? Cp’ ACD’ 
Cu Wer 


Thus, satisfaction of Ineq. (3) implies speed stability. 


7 Oo 


The Cy term used above is composed of Cp + Cas , Where Cag is 


equivalent propulsion system drag. C,. is defined as 


Gay Sees! ~ 
~ “3 3% ° 
or for constant altitude and slow speeds, 
Ga = oe 
Sav 


roportional to the negative of the thrust change 


a>) 
TS 


Cag is seen to b 


with respect to the velocity change. Inequality (3) may then be written as 


T= cet CAs _ ACd SO 


Ineq. (4) aC. 


A more sophisticated approach, based on the suppression of altitude 
disturbances by the pilot's use of the elevator, and originating with the 
basic equations of motion, would produce essentially the same inequality, 
with a small modification factor included. (See Ref. 2) The resulting 


inequality is 





Cy + Cas ACD 
Cr 


a| 








ae 7° 
Q 


However, the correction term, an ve is of relatively small magnitude 
a 


and may be ignored. From this same analysis the time constant, ae 


of the subsidence of an airspeed error is shown to be proportional to: 


Ch + Cas 
ma) De kG (SE - doy 


Typical values of the parameter I for the F-8 are given below. 
At l.l Voz, I= .22 - .0948 = - .328, showing distinct speed instability. 


pel.2V T= ,.200 - .222 = - .022, a less unstable value. 


SL! 

From an examination of the speed instability parameter, Ineq. (4), 
it will be seen that positive stability can be produced if either drag in 
trimmed conditions (Cho) or Cag is artificially increased. 

Cho could be increased, for example, by extending the speed 
brakes during an approach. This would accomplish two objectives; one, 
increase CDpo in the stability formula, and two, cause the engine to 
operate in a higher RPM range to produce the added thrust necessary, 
thus effectively causing a slight decrease in the inherent engine time 
constant. 

Remembering the definition for C age the drag coefficient corre- 
sponding to the thrust/speed variation, it will be seen that an artificial 
means of increasing thrust with decreasing airspeed accomplishes the 


desired objective. This is exactly what automatic throttle or approach 


6 





power compensator devices are intended to do. Still another method of 


increasing C_. is by the uSe of an anti-drag technique. This may be 


AS 
accomplished by having a drag producing device (such as side extending 
speed brakes) extended during the approach. Then, by using an airspeed 
controller type device to regulate the amount of drag produced (by opening 
or closing the speed brakes) in response to thrust (or anti-drag) required, 
airspeed deviations can be closely controlled. This particular method, 
within its limitations of having enough anti-drag available, theoretically 


has the advantage of a much shorter time lag than an engine response 


scheme. 





3. Artificial Speed Stability through Automatic-Throttle Techniques 

There is a choice of inputs to an automatic-throttle type controller, 
including various combinations of derivatives and integrals of airspeed 
error, pitch attitude error, angle of attack error, and normal acceleration. 
There are also many different philosophies by manufacturers concerning 
the ideal combination for best performance under all conditions. Two 
representative systems, one utilizing angle of attack and normal acceler- 
ation, and the other using airspeed error and pitch angle, were chosen for 
investigation. 

A) Angle of Attack, Normal Acceleration Inputs 

The Specialties, Inc. Automatic Power Compensator, (APC), an 
auto-throttle in use in the F-8 and F-4 aircraft, is an example of a control 
using these inputs. A block diagram, typical of this type of system is 
shown in Fig. 3. In this system, constant angle of attack is maintained 
in the landing approach using the two sensed inputs of instantaneous angle 
of attack and acceleration normal to the glide slope in g units. The normal 
g input supplies anticipation whenever the aircraft is not in a one g flight 
condition. A thrust command is calculated in the system computer depend= 
ent upon error signals and chosen gain and temperature factors. This 
thrust command then becomes an input to the engine which in turn adds 
an appropriate thrust to the airframe. This thrust, acting as the forcing 
function to the airframe modes of motion, changes the velocity of the 
aircraft. Thus the APC computer indirectly controls airframe velocity 


through the equations 





AT= Kx AK tKn Ang * AV= Kr AT 
Eq. (6) (. LV = Ky AX + Kn Ang 


where the K values are appropriate shaping, filtering and simulation 
networks. 
Or, looking at it differently, the Cae of Ineq. (4) has now been 


given a greater value ( Cas = - a 


artifically, thus making the speed stability parameter more poSitive or 
stable. 

The transfer function equations for this auto-throttle as supplied 
by Specialties, Inc., are also shown in Fig. 3. It will be noted that 
the throttle servo is simulated by a .1 second time constant first order 
lag. Angle of attack input is treated by both an integrator and lag net- 
work. Normal g input is also integrated. 

The results of the analog simulation runs using this type of control 
are reported in Section IV of this paper, Experimental Methods and 
Results . 

B) Pitch Angle, Airspeed Change Input 

Flight path angle, 3 , is given by the equation 
Ea. (7) C= O-k 

During constant airspeed flight, pitch attitude change (AG) and 
flight path angle change (4%) are proportional; but dynamically, A @—- 
occurs before A&®. Thus the means of anticipating flight path angle 


changes is available through pitch attitude changes recorded by a 


% 





vertical gyro mounted in the aircraft. This is the theory underlying the 
pitch angle, airspeed change type of automatic-throttle controller. 
Through the use of a properly designed network, these inputs make it 
possible to minimize transient airspeed changes and to eliminate steady 
state airspeed changes resulting from variations in flight path. 


The governing equation for this controller may be of the form: 


AV Kv2 AV 
Eq. (8) ke AG myn Kve AVS 
Mc =| ke (+34) 40+3% 4) 


Here the KY 2 term is inserted as a steady state error washout term. 

The lag term ( | + fsa) Simulates the airspeed sensor and filter. The 
filter attenuates the short period wind gust signals entering the system 
and acts as a smoothing influence. 

This equation was reduced to a linear open loop transfer function 
of the auto-throttle controller. The block diagram for this system, com-= 
bined with the analytically derived block diagram and transfer functions 
of the airframe and engine, are shown in Fig. 4. The stability derivatives 
for the F-8 in landing configuration were used for the airframe transfer 
functions. The derivation is shown in Appendix I. 


The system open loop transfer function was calculated as: 


f),. () 
Ue Jat (AT JAF 


Se Jaz (CC/AF \ @ ATL ATAF 


Eq. (9) 


10 





(2 uff ane - K = 
a ae AT ee AT 
ue — 62410 (1+ 2724.+18,9.4% + 235.2543 ) 


ATJAF™ 416 () 4 1.674.427. 742+ 2.08434 2,4024%) 
(£\&) = 7.014 (1+ 37.94 +176 42 
SeK“/AF 16.8 (1+5.294 +. O73 4" 
(S) — KE 
@ /Al (1+1.164) 
This equation was programmed on the digital computer to produce a 
series of Nyquist plots, each utilizing different gain constants. 
Fig. 5 shows the Nyquist plot for the equation using the same gain 
constants as were used in the analog simulation. From this figure the 
open loop gain margin is seen to be 19.2. 


CL 
The system closed loop transfer function, Sk was then derived. 


This is shown symbolically as follows: 


(At (=< 
u _ {AL JAT( AT/AF 


.(10) -—= = 
eS See” 12D) (2) eter (B or 


A digital program was devised to calculate the values of the poles 
of the closed loop transfer function for various values of auto-throttle K's. 
This procedure would have produced families of Root Loci Plots. How-= 
ever, because of lack of time and computer malfunctions, this particular 
investigation was not completed. Equation (10) would be used in an 
optimizing analysis for gain constants. 


In addition the equations of the auto-throttle were set up on the 


ll 





the analog computer. This was incorporated into the complete system 
Reine approximately optimum values of Ke , Ku, and Ku. All simu- 
lation was done in real time. These results are also reported in Sec- 
tion IV. Analog mechanization is shown in Appendix III. A composite 
block diagram showing the relative positions of the two automatic 
throttle systems in the airspeed loop is presented as Fig. 6. 

C) Auto-Throttle Controller Linked to Anti-Thrust Device 

Another possible approach to the airspeed instapility problem is 
the use of an auto-throttle controlling an anti-thrust or drag device. 
This could be Peeamolicned byflyingthe approach with speed brakes 
extended. The output of the auto-throttle would then reduce or in- 
crease the drag by closing or opening the speed brakes as required. 
Thus, instead of supplying additional thrust when the airspeed drops, 
the speed brakes would be closed a proportionate amount, reducing the 
drag and in effect accomplishing the same purpose as an increase in 
thrust. Note that the time lag of the hydraulically actuated speed brake 
system is much less than that of the engine, thus allowing a faster res-= 
——_— . The time constant of the hydraulic system is of the order of .] 
second, while that of the F-8 engine is about 1.16 seconds. This 
- approach was also investigated on the analog computer by replacing the 
engine with a first order time lag of .1 seconds (simulating the hydraulic 
system) and using the thrust command output of the auto-throttle com- 


puter to control the speed brakes. 


12 
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Obviously the limiting feature to such a system is the amount of 
drag the speed brakes are capable of supplying. Ofcourse it goes 
without saying that such a device can not be installed on an aircraft 
incapable of completing an approach while varying the positioning of 
the speed brakes. An airplane in this category is the F-8 which has a 
lower fuselage speed brake. Insufficient deck clearance is available 
when this brake is extended, thereby prohibiting its use during an 


approach. 
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4. Experimental Methods and Results 

A) Equations and Analysis 

The basic 3-degrees of freedom aircraft longitudinal equations of 
motion (in British notation) used in the analysis are shown below. The 
equations are based on the small perturbation linear theory. Lower case 
letters represent perturbation quantities. Upper case, subscripted o, 


letters represent initial steady state conditions. Body axes are used 


eee 


— 


throughout. The assumptions made are as follows: 

1. The airplane has a longitudinal plane of Symmetry. 

2. The direction of the relative wind, in steady flight, lies in 
the plane of symmetry, and in steady state, all angular velocities are 
zero. 

3. Initial flight condition is wings level. 

4. Products and squares of perturbation velocities are small and 
are ignored, and sines of all perturbation angles are approximated by 
the angle itself in radians. 

5. The longitudinal modes of motion are independent of the 
lateral modes. 

6. Structural deformations are not considered. 

Three degree of freedom non-dimensional equations of motion, 
British notation: (See Ref. 3) Figure 7 displays angular definitions and 
signs. 


X Force Equation: 
— i | — — ur ~ - a 
LH we E - [Mone th JO = Kr AT 
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Z Force _— 


Peis fS-3 ‘lie [OB | ole xi} 


Moment Equation: 


[-A"+| i * te i dalmacl oe IG Cb ma) Ie am maf 7 


where symbols are defined as follows: 


rN mM 
meus , air secs. 


d()_ ttl ) 


ip At 
‘9 = t real time 
< air secs 

uo+=« perturbation horizontal velocity, ft/sec. 


= perturbation vertical velocity, ft./sec. 

= elevator deflection, positive down, radians 
change in thrust required, lbs. 

ze XK , angle of attack, radians 


= pitch angle, radians 


I 


-@ gle Be * 


ne — a2 oe 

b m Le2 , inertia Coefficient 
Joy = Moment of inertia, slug = ft 
mM = mass, slugs 

4, = tail length, ft. 
A an — = relative densit 

1 = 
PS 2 e Y 

é - mean wing chord, fi. 

S = reference Wing area, square it. 
P= Force parallel to x - axis, lbs. 
a = Force parallel to Z — axis, lbs. 
MM. = Moment about lateral axis, lbs. = ft. 
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Kew 


Xr 


Mu 


distance from thrust line to CG, ft. 


ve 
Oj 
1/2 (j Cus non-dimensional Force Stability derivative 


SS 


&) 


and similar identities obtained hy permuting 
xy 2 and u,& independently. 


ax: 
1/2 al =) non-dimensional stability derivative and 


2 (Fs) 


WZ fas 


enn 


ove 


Fes) 
W2Ct Gu je 
*3) 


x 
WV (Zeus) 
>o 


ne eos eS) 
it y(ucs) 
Ue 


Xr 3 


Similar identities obtained by permuting x, 


a 


non-dimensional stability derivative and 
similar identify obtained by permuting x 


and 3. , n and T 


non-dimensional moment stability derivative 
and similar identities obtained by permuting 
uand wr 


non-dimensional stability derivative and 
similar identities obtained by permuting x, 


and ee 


non-dimensional stability derivative 


non-dimensional stability derivative h = 
distance from thrust line to C.G. position 


The necessary stability derivatives were evaluated by making use 


of a longitudinal dynamic stability Fortran program, designated LONGSTAB. 


The program was compiled while undertaking the course in dynamic air- 


craft stability given at the Postgraduate School. This program computed 


the necessary British stability derivatives and then used these derivatives 


in the aerodynamic equations of motion. The resulting stability quartic 


was then solved, by the program, for the phugoid and short period 


complex roots. Finally, the periods and times to damp to 1/2 amplitude 
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were computed for both modes. Necessary aircraft parameters uSed by 
the program were taken from F-8 data supplied by Ling-Temco-Vought, Inc. 
or calculated. The program-computed stability-derivatives used in the 
investigation are shown in Table I, as are the other aircraft parameters. 
All program-computed derivatives, after necessary conversions were 
found to be in close agreement with those listed in Ref. oa » The computer 
program is explained in Appendix II. A sample printout is shown in 
Table II. 

From the program printout, the phugoid period is 33.5 seconds and 
that of the short period mode is 6.06 seconds. Times to damp to 1/2 
amplitude are 29.9 and 1.73 seconds respectively. Undamped natural 
frequency , co Bmagoide as calculated as .189 radians/sec., 5S asosl 23, 
Short period G)_ is calculated as 1.19 rad.Jsec., J as. .360. 

Ref. 6 lists a phugoid W, of 0.188 rad./sec., anda Yof .12. 
Short period results found were listed as CJ = 1.19 rad./sec., § = 
0.35. The results given by LONGSTAB and those given by Ref. 6 are in 
excellent agreement. 

Thus, as previously remarked, it is seen that the phugoid is very 
lightly damped. 

lref. 6 reports on the results of an investigation conducted along 
Similar lines as the subject of this paper. The two investigations were 
performed concurrently and independently of each other. After com- 


pletion, the results of the author’s analysis were verified insofar as 
possible with those given in Ref. 6. 








These non-dimensional equations of motion were dimensionalizecd 
for the analog simuiation. The engine time lag and auto-throttie blocks 
were next added. Real time stick and throttle simulators were then in- 
troduced into the system, together with the F-8 linearized longitudinal 
control system dynamics supplied by Ling-Temco-Vought. All analog 
runs were conducted in rea] time. The analog computers used were 
Donner 3100 and 3400 machines. An eicht channel Brush recorder was 
used to record the time histories. 

The followins types of computer runs were tnen conducted: 

(a) Transient analysis runs using step inputs of elevator deflection, 
horizontal velocity gusts, and pitch attitude: first with basic airframe 
alone and then with the two basic types of auto -ihrottle controllers in- 
corporated. Finally sample runs were made with the speed brake device 
replacing the engine. 

(b) Optimizing runs using triangular inputs of elevator deflection, 
to determine the optimum gain tor this type of forcing function. 

(c) Sinusoidal gust inputs. 

(d) Real time pilot controlled runs using a simulated glide slope, 
with and without the aid of the auto-tnrottle. Horizontal turbulence was 
also added to determine the controller efiectiveness. 

Analog mechanization diagrams are shown in Appendix III. Poten- 
tlometer settings are shown in Table 1 to Appendix III. Time histories for 


these runs are displayed in Figs. 19 to 35. 





B) Evaluations of Experimental Analog Runs 


The systems referred to are as follows: 


System 1, Inputs of angle of attack, and normal acceleration 





in g units. 
System Equation: AT =—! Kx 4% _ KavAV Ky om | 
It IAL [RSA | +A. .. 


System l-a gain constants; 


Kx = 1970 lbs./degree 
Kav = 280 lbs./knot 
K, = 183 lbs./degree — sec. 


System l1-b gain constants; 
Kx = 3420 lbs./degree 
Ke =. (2 les knee 
K, = 190 lbs./degree - sec. 
System 2, Inputs of pitch angle, airspeed error 


-“&_ [| Ku+ ku 
System Equation: AT = KoA el aa A. | 


Gain constants; 


Ke = 318 lbs./degree 
Ke = 360 lbs knot 
Kjc = 22.3 lbs./knot - sec. 


1) Step input of 1.32° up elevator deflection 

Fig. 10 is the analog time history response toa 1 ade step of up 
elevator deflection, using the System l-a auto-throttle. Airspeed dis- 
plays a slight oscillatory tendency, reaching an equilibrium speed of 
+ .6% of ee (.8 kts. above V) » Thrust is applied smoothly in response 
to demand. 

Fig. 11 shows the response of System l-b. This system had the 
gain factors optimized for a step input of horizontal gust velocity. Note 


the high airspeed overshoot due to step elevator input. 
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nirspeed is much 


Fig. 12 is the same situetion usiry System 7 
more highly damped as compared to Fix. Eouilibrium airspeed is 
-1% of Ve (1.39 kts. below V,). Thrust is also applied smoothly. 
2) 5 knots of horizontal gust input 
ste response to this 
tamping ratio of 


oversnoct is 61% of 


‘ 
@ 


par 


Fig. 13 displays the recora of S 
ess than optimum for 


input. The system is oscillatory with ar approximate c 
? he Q 


, which has been 


5 = .15,@) = .677 rad/sec. 1 
These values of wain co 
now 98, Od, is 


Ps 
~~ 


) 


ra Lex 
Fine! value of 


initial error. 
this type of input disturbance. 


Fig. 14 displays the response o% s 
nat oscillatory. 


optimized to this input disturbance. 
.73 rad/sec. Thrust response is somew 
airspeed is +.2% of - : 
Fig. 15 is the record of Srstem “©. 1% shows @ damped response 
with final airspeed of about =% .0%. ° Loriet changes are less than 
that required by System |. 
3) 4° step input of pitch aniie 
Fig. 16 shows the response of S$, S*er *o this input. Initial 
velocity transient is +.&% a or i 12 ete. steady State error is +.4% Vo: 
Thrust application is slightly o¢c!)lacor 
l-r yvesponse. Initial velocity transient is 
i ny 
G 
initial velocity transient 


Fig. 17 shows System | 
steady state error in 193 


.6%V . Steady state error is also .%% 
oO 
Fig. 18 displays the resport se 0% 5°"s5* 
Airspeed damps *o +.4 


amounts to +1.6% V 





seconds, reaching less than 1% error in 9 seconds. 

Fig. 19 shows the response of the basic airframe to this same 
input disturbance, without the ald of the auto-throttle. Airspeed 
oscillation is +1% Vos with a steady state value of about .5% wee 

4) Optimizing runs using triangular elevator input 

Fig. 20 shows System l=-a under the influence of a .02 cps 
triangular elevator deflection (maximum amplitude 1.32°) forcing 
function. Airspeed error is kept within a rance of +2% to -1% of V 0° 
Thrust varies within the range of +1600 to -2400 lbs. about the trim 
setting. Thrust application follows elevator movement closely. 

Fig. 21 shows the reaction of System l-b. Note how airspeed 
fluctuates + 4% we (+ 5.56 kts.). Thrust has very noticeable oscillatory 
tendencies and gyrates within rather wide limits. 

Fig. 22 shows System Z2 under the influence of this same forcing 
function. Airspeed is seen to vary + 1.7/7 kts. ort i.c%ofV,. Thrust 
varies approximately 2400 los. about trim thrust. Application of thrust, 
however, is smoother and displays less hunting characteristics than 
shown by System l-a. 

Figs. 23, 24 and 25 show the response of these same systems to 
@ triangular elevator forcing function of the same magnitude as previously, 
but of .06 cps frequency. 

Fig. 26 shows the reaction of the airplane to this same forciny 


function without the benefit of the autc-throttle. Note airspeed deviation 


is less than that obtained with System 1-b engaged. 





S) Human pilot controlled flight on simulated glide slope 

For this series of runs a pilot fiew the simulated eirplane on the 
glide slope, using the auto-throttle. Glice slope reference was obtained 
from a zero centered panel meter. The zero voltage level was considered 
the correct glide slope. This is shown on Channel 1 of Plgs. 27, 28 and 
29. The pilot was then told to manipulate the stick to fly the aircraft to 
a new A higher glide slope. This corresponded to a+ 14 volt reading on 
channel 1, or 14 volts on his reterence meter. This simuiated a low 
"meatball" on a mirror approach and the flight path correction necessary. 
This new flight path angle was held momentarily. The pilot was then told 
to fly the airplane down to a = Ae glide slope, aiter which he was to 
return to the normal flight path (0 volts). 

The response of the aircrait, and speed control is shown in Figs. 
27, 28 and 29 for Systems l-a, l-b and 2 respectively. 

It will be noted that Systems i-a ana 2 mainiein sirspeed close to 


within 1% of V_, whereas System I-b has much areater airspeed fluctua- 


an 


° 


tions. Thrust application is again much more oscillatory in System l=b. 


6) Human pilot controlled flisnt on alide slope with sinusoidal 
input of 5 knots horizontal velocit: Just 
Figs. 30 and 31 are records of System 1 runs utliizing a pllot 


« 


controlled stick with a forcing function sinusoidal input of 5 knots of 

horizontal velocity, frequency .05 cys. The pilot was again attempting 
O 

to fly the glide slope 4 up and down. 


Fig. 32 is the record of the response to this same forcing function 


under the influence of the pilot and the Systern 2 auto-throttle. 





Fig. 33 displays the record of the pilot attempting to maintain 
original glide slope while the aircratt was being influenced by a sinu- 
soidal u turbulence forcing function, without the benefit of the auto- 
throttle. 

7) Speed brake anti-thrust system 

Fig. 34 shows the history of the speea brake anti-thrust system 


in response to a 5 knot horizontal velocity step input. Mechanically, 


this thrust command was used to drive @ .! second time constant antl- 


drag device, simulating the movement of the speed brakes. It will be 


4 


>, 


noted that in comparison with Fig. 15, (standard System 2 setup) response 


+ 


and damping is faster. It required 9 seconds for the transient to steady 
within 1% of es in the case of Fig. 34, while it required 11 seconds for 
System 2. 

Fig. 35 shows the response of the anti-thrust system to the same 
forcing function, but utilizing different gain factors. Response is more 
oscillatory. 

C) Summary of Analog Resu!ts 

From the foregoing evaluations it can be seen that the response of 
the different systems investigated is dependent upon the sensed variables 
and the applied forcing functions. System l-a, optimized to an elevator 
deflection input, is marginal to TRIBES enor in performance when 
subjected to a horizontal gust input. On the other hana, system l-b, 


using the same sensed variables, but with gains optimized to horizontal 





velocity input, does not perform well wien subtected to elevator 
deflection inputs. System 2, utilizing different sensed inputs, also 
displays varying reactions to each oi the forcing functions. 

In general, it can be seen that System 2 has the capability of a 
smoother application of thrust in response to an airspeed error. System 1 
has the capability of a more highly damped response, but causes a more 
erratic application of thrust. Further investioation as to the optimum 
values of gain constants is required to accomplish the desired task more 
efficiently. 

It is further recommended that additional sensed variables or 
integrals of such terms be added to the analog set-up to help optimize 
system response. A fruitful area miaht be to incorporate a normal g 
input to System 2 to help the damping characteristics. It is anticipated 
that further investigations along this line will be conducted at the U. S. 
Naval Postgraduate School. These investigations will attempt to utilize 
the simplified equations obtained by the elimination of the short period 


mode, to develop a more tractable analysis. 
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oe Conclusions 

The two major approach power compensator systems investigated 
appear to have the capability of performing the aesired task. The major 
difference are degree of anticipation and erraticism ot thrust application. 
With proper choice of gains and input variables, it seems feasible that a 
composite system can be desianed to provide the advantages of both 
individual systems. 

The speed brake or anti=-thrust system has the capability of pro- 
viding faster response, but is limited by dray area available and the 
physical feasibility of approaching with speed brakes extended. 

The preceding investigation has firmly convinced the author that 
an auto-throttle mechanism is an important contribution to carrier aviation 
safety. An aircraft that is difficult to tly in tie approach because of 
inherent speed instability problems, can hawe its whole character changed 
by an auto-throttle. By removing the undesired instabliity an auto-throttl 
relieves the pilot of concentrating on airspeed control, this allowing more 
precision in glide slope and line~up control, Tris advantage by itseli, 
should tend to reduce the frequency of ramp accidents. 

In addition, when an automatic landing system is considered, the 
pilot-controlled airspeed loop musi be replaced by an automatic system 
to provide consistency and predictability of response. 

Still another advantage of the automatic airspeed control system is 
the reduction of minimum approach speed helow that recommended by the 


pilots when airspeed is controlled by manual throttle manipulation. In 





addition the automatic device can he desizres to provide the anticipation 
and instant response that can only be chtaines from nricn actual pilot 


experience. Thus it is seen that the aute-throttle can act as the great 


equalizer in carrier approach trainin. 
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Computed F8 Stability Derivatives and Associated Aircraft Parameters 
British Non=Dimensiona! Notation 
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INPUT PARAMETERS 
TAIL AREA= 95-nCOC WING ARe b= 275.C000 CL ALEA TAIL= 2.293CCC ETA TAIL= 1.0C000U RCLY AREA= be. 20 
BODY LENGTH=52.80C00 RHC= .0C22780 GEE=32017 080) wEIGHT= 2-Co0.cGCC Ces D EPSILCN/D ALFA= 477200 
INITIAL PITCH ANGLE= €.10C000 GEGREES THRUST ee pie heheee =.437€0 FEET TERUST CCEF CT= -195CCC 
CMALFA= --28000C eave y= GéECCO-CCO SLUG FT SQUARED 
STABILITY OERIVATIVES 
FORWARD VELOCITY= 234.00C T CARAT= 302712735 G= 65-1049 
C SUBL= ~90000C C SUPO= ~26300C 
XU= -.2630 XW ~ 16996 XQ= ©COOO XTRETA= --4S5000 ZL= -.9C000 Zh* -1-5915 TOTAL ZQ= -.5878 
2ZTHETA= --0640 WU= .C06052. MW= --15896 MBAR wOGTs <a bh22 C= ~-f£ti9 MUCNE= SHHE6 I SUB BETA= ~7Cel 
A= 1.00000 8= 2e71 ech C= 14.0R4677 C= R.CZGz2e SEs 5.Caush 
THE VALUE OF ROUTES DESCRIMINANT = 69.47 
THE AIRCRAFT IS STABLE 
RUN 3 CeGe 245 V2139 KTS,234 FPS 3/3/62 LEVEL FLITE 
THE COEFFICIENTS CF THE GIVEN PCLYNCMIAL ARE 
REAL PART 
-0000000 2.776C450 14 .0467662 3.0202908 5.0845293 
IMAGINARY PART 
-0000000 -ooocooc ~cCccooo ~occccce eccccccc 
’ TABLE OF ROCTS 
ROOT REAL PART IMAG PART PERICC TIME TC CAMP TO FALF AMPL 
NUMBER OF RcOoT OF RCCT 
1 -~.C75920 - 614444 2a, £1208 SECS 29.9173 SECS 
2 -.C€75930 2 614bub az cy7e SECS 29.9173 SECS 
3 - 1.312092 -3.397553 6.C6C¥% SECS 1.7313 SECS 
4 -1.312052 32397553 6-C6C? SECS 1.7313 SECS 
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Where: 
AT = Thrust Command 
AX'= Angle of attack err. 
M%, = Reference. Angle o + 
QO = True Ang le of atts 
Anz zs ch enge ot normal OW 
Vo > Reference Airs p eed 
Ky Angle o€ attack gels 
Ky S peed Change Gain 
Kx Integ ral of angle of all Lin actor 
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APPENDIX I 
Derivation of Transfer Functions 
The dimensional derivatives shown in Table I were usec to calcuiate 
the analytic transfer functions. Calculated values of these cerivatives 
using F8 data, are shown in Table II. 
The longitudinal Equations of Motion written in Matrix Laplace Form 


are (neglecting negligible terms); 


a-Ku “Ky 79 (4) O Kar\ fde 
oe 4-Z7ur mans OF (a) = Ze O 
“Mu “4M Mus) aM G(4) Mse Mat AT 


ve . 
The Transfer function ( fe Airframe then can be written as 
Q 


(in determinant form): 

oO AG amg 

Bie A-Kuw -aUp 
Mse -(aMietMus) 4° Me a 







D, 7) 
Where D, = determinant of homogeneous equations 
Expanded D, = A 4* +PaesC 2cueemey F 
Where A = 1 
B= -(Me tXu +Zur *UoMur) 
c= M¢tu, Us Mist Xu (Me t2urtU, » Gs aia 
D= XE (Mg Zu Vou)" MuUoXut pK B,79(MicBi Mu) 
25 9(Muzu =f tee) 
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Expanded Nee = Bee An + Gore: + Du 
Where Bu= Zee Xur 
“Fre (gMir+Mag Xw) +Mse (Uo Kur 74) 
Du = G(Mse Zur“ Mw Zse) 


Substituting values for these cerivatives 


a _ waif 424 4504 785.2] _ Ne 
Sehvtrame 22+ 86547 +11.52947 4.69644 416 DI 


By a similar procedure, 


(£-\ _ —2.50[ 42+. 21554 +.00568] 
Se JAipframe D, 


(Positive elevator deflection is down) 


The determinant for calculating (< an »@ written as 


= 
Xat "Xu TS 
© A= ow ~AV 


Mar -(@Mus +Mu-) 4*%Mq 
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Expanded: 
fe — Arad +Bra?+Cra Dr 
AT )irtrame. Dd, 
Where 
Ar= Xar 


By= “Xar| Zu-+Me|+ Uo Mur 
Cp = Xarl 2M, “UoMw] u Mat [Ue Xu 9] 
Dr = Meat q ae 


After substituting in values: 
(25 —- 200145 Ca? + 80994 ih S94 1.0430] 
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APPENDIX —_L 
TABLE I-A 
LONGITUDINAL DIMENSIONAL STABILITY DERIVATIVE 
PARAMETERS 
(STABILITY AXIS SYSTEM) 


a ee Re ee A GR 


IN TERMS OF BASIC STABILITY DERIVATIVES IN TERMS OF NON- 
—_——_—_——— a DIMENSIONAL 
QUANTITY DIMENSIONAL NON -DIMENSTONAL STABILITY DERIVATIVE 
== PARAMETERS 
DEFINITIONS UNIT 
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TALLE J-B 


LONGITUDINAL NON-DIMENSIONAL STABILITY DERIVATIVES 


(STABILITY AXIS SYSTEM) 





BASIC NON-DIMENSIONAL STABILITY DERIVATIVES NON-DIMENSIONAL 
= 
TOTAL AIRFRAME THEORETICAL HORIZONTAL  . STABILITY DERIVATIVE 
DEFINITIONS UNIT TAIL CONTRIBUTION PARAMETERS 
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App end IT fable IE 
Dimensional STabihty. _Derivetives 
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APPENDIX II 
DYNAMIC STABILITY FORTRAN PROGRAM 

Longitudinal - Program LONGSTAB 

This program, once entered with required airplane parameters , 
computes the longitudinal stick fixed and/or stick free (when pertinent) 
stability derivatives, and the corresponding stability equation coeffi- 
cients. It then solves these equations for roots, periods, and times to 
damp to 1/2 amplitude for both the phugoid and short period modes. In 
addition the value of Routh's discriminant is determined. The sequence 
of computations can be repeated for as many different flight regimes as 
desired. 

There are a total of five data cards required for the stick-fixed 
only solutions, with an additional two cards required for a stick-free 
Paintton. Thus for a combination stick-fixed and free solution, there 
will be a total of seven data cards required for each flight regime 
desired. 

The initial data card controls the types of solutions desired 
(stick-fixed, free, or both), and the number of consecutive runs to be 
made. A list of program symbols with their meanings appear in Table I 
to this Appendix. 

The second card contains an alpha-numeric run identification or 
an arbitrary title. 

The third card is the first of the general input data cards. This 


card should have entered on it: 


eZ 








S, (tail area), S_ (wing area), Chet Ve tail’ Sp (body area) , 1, 
(body length) , g , g(acceleration). 
The parameters are entered eight to a card, each being allowed 10 
columns. There are 80 columns on each card. Each parameter may 
be entered anywhere in its assigned ten columns, but must have a 
decimal point included in it. This same format is followed for each 
of the succeeding cards. 

The fourth card should have the following aircraft parameters 


entered on it: 


ne 
3k 


c.g. to thrust line, ft.); C. (coefficient of thrust): C 


W (aircraft weight, lbs.); ; CE in degrees; h, (distance from 


Z 
ye Iyy (slug ft .): 
V (forward velocity ft./sec.). 

The fifth card should have the following: 


Cp: Cre! CITA: C * c (wing chord): ” (eyes tora. ec. tail sft. ): 


L wing 
Xc.g. (x distance of c.g., ft.); XAC (x distance of wing a.c. in ft.). 

Cards 2 through 5 are repeated for each different run, stick fixed. 
If it is desired to compute a stick free solution simultaneously, then 
two additional data cards (6 and 7) are required for each run. 


Card 6 contains the following parameters: 


Cr (tail chord, ft.); Ke (radius of gyration, squared) Xe ; By: Bo: 


gi Bos 


re? 7: 


Card 7 contains the following parameters: 


e e e o 2 
My + My im, (mass elevator, slugs); S. (area elevator, ft:). 


3 





Thus for one aircraft flight regime (say Sea Level, Mach 0.6), 
for both stick fixed and free, there would be required data cards 
1 through 7. For each succeeding run, cards 2 through 7 would have 
to be repeated with the new parameters entered. Card one would have 
to be suitably prepared to reflect the total number of runs and mode of 
runs desired. See Fig. 1 this Appendix for sample data cards. Note 
the first card is of format 13 for the number of runs, and I! for mode of 
operation. All other data cards are of 8F10.0 format (each card holds 
8 fields of 10 columns, each, one data number to be entered per field, 
in floating point format: i.e., with a decimal point somewhere in the 
number.) 

The program prints out the calculated stability derivatives as 


well as the computed solutions to equations. 
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TABLE I 
PROGRAM SYMBOLS AND MEANINGS 


Program LONGSTAB 


Program 

mnenonic Meaning 

NOS Number of runs to be made. I3 format. card 1 

MODE Decision stick fixed (1), or stick free (2). Il. card 1 
ST S;, tail area. card 3. 

SW Sys wing area. card 3. 

CLALFAT C , lail lift curve slope. card 3. 


L&T 
ETATAIL N, » Tail efficiency, C,/Cy. card 3. 


SB Sp: Body area card 3. 

BL 1a: Body length Cdrdsc. 

RHO © , Atmospheric density, card 3. 

G g, Acceleration due to gravity. card 3. 
W w, aircraft weight. card 4. 


EPSALFA 9£ card 4. 
aK 


THETAI 6, , initial angle of pitch, degrees. card 4. 
YH h, distance from c.g. to thrust line. card 4. 
CT Cre coefficient of thrust. card 4. 
CMALFA C 2 per radian. card 4. 
MX 
YI ly! moment of inertia. card 4. 
V V, forward velocity, ft./sec. card 4. 
CSUBD Ch: coefficient of drag. card 5. 
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CLALFA 


EPIA 


CSUBL 


TL 


XCG 


XAC 


CTAIL 


XKE 


UE 


Bl 


B2 


B2PRIME 


ZETA 


ZETADOT 


AMETA 


AMETADT 


XME 


SE 


Shi , Winglift curve slope. card 5. 


CITA, card 5. 
Cy wing’ card 5. 
c, wing chord. card 5. 
L, C.g. tO a.c. tail. Sear... 
ee , X distance of c.g. card 5. 
Xo. X distance of wing a.c. card 5. 
Cy, Tail chord.  cards6- 
Re radius of gyration squared. card 6. 
Fre relative elevator density. card 6. 
B., C card 6. 
Lo her 
Bas Chee Cara b- 
Bo, 1/2 Cup , card 6. 
a "Ceara 6. 
ar! card 6. 
Wieeescards7 . 
ne 
Mm: , card 7. 
aE 


moe mass elevator. card 7. 


Se: area elevator. card 7. 
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APPENDIX IT! 
SYSTEM MECHANIZ ATION FOR THE ANALOG COMPUTER 

Fig. 1 to this Appendix displays the analog schematic for the 
equations of motion of the airframe. Fig. 2 includes the schematic 
for the System I auto=throttle, aircraft engine, and stick system 
dynamics simulations. Fig. 3 shows the System 2 ute aera 
schematic. 

Table I lists the potentiometer settings for all systems. 
Table II lists the scale factors used in computing these settings. 


Table III describes the calculations used in arriving at these values. 
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TABLE tI 
SCALE FACTORS USED IN ANALOG SIMULATION 
Scale Factor Value 
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